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ABSTRACT 


An  analysis  is  presented  of  two-dimensional,  isoenergetic,  com¬ 
pressible  mixing  of  a  jet  with  a  fluid  at  rest  for  both  laminar  and  tur¬ 
bulent  mixing.  The  analysis  is  primarily  concerned  with  the  development 
of  theoretical  expressions  for  the  mixing  similarity  parameters.  A  gen¬ 
eral  momentum  equation  is  derived,  which  relates  the  similarity  param¬ 
eter  to  the  mixing  length  and  viscous  shear  stress.  By  using  this  general 
equation  and  Newton's  viscous  shear  relation,  a  complete  theoretical 
solution  for  the  laminar  mixing  similarity  parameter  is  derived,  which 
does  not  involve  a  reference  perturbation  velocity  factor.  Prandtl's  mix¬ 
ing  length  theory  for  the  apparent  viscous  shear  relation  is  used  to  obtain 
the  theoretical  similarity  parameter  for  turbulent  mixing.  Based  on  these 
similarity  parameters,  equations  for  the  width  of  the  corresponding  mix¬ 
ing  regions  are  derived.  These  results  and  the  approximate  theoretical 
velocity  profile  equation  developed  by  Pai,  Nash,  and  Korst  represent  a 
closed-form  theoretical  approximation  of  the  type  of  mixing  considered. 
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SECTION  I 
INTRODUCTION 


The  type  of  mixing  considered  is  two-dimensional,  compressible, 
isobaric,  and  isoenergetic  without  an  initial  boundary  layer.  The 
object  is  to  define  theoretically,  in  closed  form,  the  following  character¬ 
istics  for  both  laminar  and  turbulent  mixing: 

a.  The  velocity  profile 

b.  The  width  of  the  mixing  region 

For  laminar  mixing,  both  of  these  characteristics  were  numerically 
obtained  for  specific  cases  by  Chapman  (Ref.  1).  An  approximate  velocity 
profile  for  laminar  mixing  was  obtained  by  Nash  (Ref.  2),  who  used  the 
simplified  equation  of  motion  of  the  heat  conduction  form,  derived  by  Pai 
(Ref.  3).  However,  this  equation  contains  an  unknown  reference  perturba¬ 
tion  velocity  factor  for  which  Pai  assumes  a  value  of  unity  and  Nash  assumes 
a  value  of  2.  0.  Nash  verifies  his  assumption  by  comparing  with  the  numeri¬ 
cal  values  obtained  by  Chapman.  For  turbulent  mixing,  Korst  (Ref.  4)  also 
reduced  the  equation  of  motion  by  the  method  of  small  perturbations  to  the 
heat  conduction  equation  form.  Thus,  the  velocity  profile  for  both  laminar 
and  turbulent  mixing  is  approximated  by  the  equation 

</>  =  -J  [!  +  erfTj] 

where  the  error  function 


The  difference  between  laminar  and  turbulent  mixing  is  reflected  in  the 
relationship  for  the  similarity  parameter,  a  .  For  turbulent  mixing, 
experiment  shows  a  to  be  a  function  only  of  the  free -stream  Crocco  number 
since  the  mixing  zone  width  varies  linearly  with  the  mixing  length,  X  .  For 
laminar  mixing.  Page  and  Dixon  (Ref.  5)  developed  the  following  expression 
for  a  (which  can  also  be  obtained  by  inference  from  Nash's  work): 


where  a  is  the  unknown  reference  perturbation  velocity  factor.  Thus,  in 
general,  the  similarity  parameter,  a ,  is  a  function  of  both  the  free-stream 
conditions  and  the  mixing  length,  X  . 

The  following  is  a  derivation,  which  yields  a  complete  theoretical  rela¬ 
tion  for  the  laminar  mixing  similarity  parameter  and  a  semi-theoretical 
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relation  for  the  turbulent  similarity  parameter  based  on  Prandtl's  mix 
ing  length  theory. 


SECTION  II 
BA5IC  EQUATIONS 


The  purpose  of  this  analysis  is  to  relate  the  mixing  similarity  param¬ 
eter,  a ,  to  the  mixing  length,  X  ,  and  the  viscous  shear  stress,  r,  in  a 
manner  that  is  applicable  to  either  laminar  or  turbulent  mixing.  Such  a 
relationship  can  be  derived  by  applying  the  conservation  of  momentum  con¬ 
dition  to  the  control  volume  shown  in  the  following  figure  of  a  general  mix¬ 
ing  zone. 


1.  Over  the  mixing  zone  length,  X  ,  considered,  the  average 
angle  between  the  dividing  streamline  and  the  mixing  zone 
centerline  is  small  (<20°). 

2.  The  mixing  is  two-dimensional,  isobaric,  and  isoenergetic 
with  the  non-dimensional  velocity  profile  independent  of  X  . 


The  momentum  equation  for  this  control  volume  is 

x  yd 

/  TDdX=  J  pii  dY 

n  *■  V  >f 


(1) 


By  the  Perfect  Gas  law 


P  = 


p 


RT  RTt{l-CHM^) 
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where 


let 


where 


*  - 


D 


U 


00 


=  1  - 


Tt 


<7  =  similarity  parameter  and  is  a  function  of  Ca  and  X  . 

eo 


Substituting  into  Eq.  {1)  yields 


Let 


(2) 


The  integral,  I3,  has  been  numerically  evaluated  in  Refs.  6  and  7  based  on 
the  non-dimensional  approximate  velocity  profile  discussed  in  Section  I. 

As  a  consequence  of  assumption  2,  I3  is  independent  of  X  .  Therefore,  the 
derivative  of  Eq.  (2)  with  respect  to  X  is 


rD  = 


<3) 


Equation  (3)  is  the  desired  general  relation  between  the  viscous  shear 
stress,  rD  ,  the  mixing  length,  X,  and  the  similarity  parameter,  a.  Equa¬ 
tion  (3)  also  embodies  the  conservation  of  the  X-direction  momentum  and 
mass  flow  in  the  transverse  Y  direction  since  the  location  of  the  dividing 
streamline  is  determined  by  applying  these  conditions. 


2.1  LAMINAR  MIXING 


The  similarity  parameter,  a,  for  laminar  mixing  can  be  derived' 
from  Eq.  (3)  in  the  following  manner.  For  laminar  mixing,  the  viscous 
shear  stress  is  determined  by  Newton's  relation  as  follows 


rD  ~ 
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Since 


-($)(&)-■-(£)(£) 


therefore 


J2L_  _ 

dY 


rD  =  u« 


Substituting  Eq.  (4)  into  Eq.  (3)  yields 


-o  O.G)  -  [fr  <hd] 


Separating  the  variables  and  integrating  yields 


T [*  '«,]  (S' 


Solving  for  a  gives 

The  width  (b)  of  a  laminar  mixing  zone  can  be  determined  from  the  definition 
of  7j  and  Eq.  (5)  as  follows  since 


•  *  b  -  2tjm  -£~ 

Substituting  Eq.  <5)  for  o  yields 


b  =  2 


~ 2  BT-  *■  (4£)„  »" 

p  u»  (i,) 

*?D  J 


Thus,  for  given  free-streapi  conditions,  the  width  of  a  laminar  mixing 
zone  varies  linearly  with  X  ^ .  The  value  of  7m  is  determined  when  the 
limits  of  the  mixing  zone  are  defined.  A  mixing  zone  is  usually  defined 
to  be  that  region  which  exists  between  7M  =  ±1.503  ,  which  corresponds  to 
velocity  ratios,  0  ,  of  0.  9845  and  0.0154  based  on  the  approximate  velocity 
profile  equation.  This  definition  is  based  on  Tollmien's  results  in  Ref.  7. 

From  Eq.  (5),  the  7  ordinate  for  laminar  mixing  is 


7  =  ° 


Y  pu»U 

-X  =  Y  - 

X  / 

_2RTt  ( 


(—) 
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since 


Substituting  into  Eq,  (7)  yields 


In  Ref.  1  velocity  profiles  for  laminar  mixing  are  presented  as  a  function 
of  v  (  u°°—  )  and  are  shown  for  various  free-stream  Mach  numbers  in  the 

'  V 


range  from  0  to  5.  0.  The  velocity  profiles  for  free-stream  Mach  numbers 
of  0  and  5.  0  have  been  replotted  in  Fig.  2  as  a  function  of  jj  (Eq.  8)  using 

the  compressible  value  of  (1,)^  from  Ref.  8  and  evaluated  from  the 

approximate  mixing  velocity  profile  also  presented  in  Fig.  2. 


t 


Fig.  2  Correlation  of  Laminar  Mixing  Velocity  Profiles 
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The  parameters  were  evaluated  at  the  dividing  streamline  location  shown 
in  Ref.  1  to  correspond  to  <£D  =  0.61  for  n  =  0.76  .  Figure  2  shows  that 
the  velocity  profiles  for  laminar  mixing  are  correlated  very  well  by  the 
V  variable  defined  by  Eq.  (8).  The  deviation  of  the  exact  velocity  profile 
for  M  =  5.0  from  the  approximate  theoretical  profile  is  due  to  the  fact 
that,  for  laminar  mixing,  the  exact  velocity  profile  is  a  weak  function  of 
the  free -stream  Mach  number.  This  is  also  indicated  by  the  invariance 
of  the  dividing  streamline  velocity  ratio  with  free-stream  Mach  number. 


In  Ref.  5  for  laminar  mixing,  the  following  relation  for  a  is 
presented 


(9) 


where  a  is  a  reference  perturbation  velocity  factor.  A  theoretical  rela¬ 
tion  for  a  can  be  obtained  from  Eq.  (5)  in  the  following  manner.  Eq.  (5) 
can  be  written  as  follows: 


(10) 


Equating  Eqs.  (9)  and  (10)  and  solving  for  a  yields 


The  theoretical  variation  of  a  with  free-stream  Mach  number  is  presented 
in  Fig.  3.  The  parameter,  a  ,  is  shown  in  Fig.  3  to  have  a  value  of  2.  112 
for  Ca^  =  0  ,  which  agrees  with  the  value  of  2.  0  assumed  by  Nash  (Ref.  2). 
However,  o  increases  in  value  as  the  free-stream  Mach  number  increases. 


Fig.  3  Theoretical  Variation  of  the  "Reference  Perturbation  Velocity  Factor'* 
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2.2  TURBULENT  MIXING 


The  similarity  parameter,  a  .for  turbulent  mixing  can  be  derived  from 
Eq.  (3)  in  the  following  manner.  For  turbulent  mixing,  the  apparent 
viscous  shear  stress  can  be  defined  by  Prandtl' s  mixing  length  theory 


rD  -  Pn  & 


°r  Pn  <J1  ’’  (±£\  U2) 

r°  '  P  XJi/d 

Since  the  mixing  length,  l,  is  linearly  related  to  the  effective  vortex  size, 
then  the  following  relation  for  £  is  consistent  with  assumption  (2). 


&  =  kb 


where 


k  -  non-dimensional  unknown  constant 


By  definition 

gb 

%  2  x 

or 

b 

■■  l  -  2  TJy  k 

Since  S.  varies  linearly  with  the  ratio,  X/g  (not  X  as  assumed  by 
Tollmien  and  others),  the  effect  of  compressibility  is  taken  into  account 
and  the  constant,  k  ,  will  be  independent  of  the  free -stream  Mach  number. 
Substituting  into  Eq.  (12)  yields 

rD  =  ^  k 2  Pd  u«  (^)d 


Pd  htd  RTt(l-C4^D2) 

r  _  M" 

D  RT‘(1“C^^d)  \dl?/D 

Substituting  Eq.  (13)  into  Eq.  (3)  yields 


Separating  the  variables  and  integrating  yields 


(13) 


(14) 
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Tollmien  (Ref.  7)  experimentally  determined  a  value  of  12  for  a  from  low 
subsonic,  Ca^  »  0,  turbulent  mixing  experiments.  Based  on  this  result 
and  using  the  approximate  velocity  profile  and  the  theoretical  dividing 
streamline  location  from  Ref.  8,  the  value  of  the  constant,  k  ,  in  Eq.  (14)  is 
0.06895  for  =  1.503.  Substituting  these  values  into  Eq.  (14)  yields 


«^d(1  -  c.-  *d) 

(  0.04235)  \ 

\*v)a 

The  equation  for  n  relative  to  the  incompressible  value,  Ca  =  0,  is 

oc 


(15) 


a 


o 


a 


(16) 


Fig.  4  Compressibility  Effect  on  the  Similarity  Parameter  for  Turbulent  Mixing 


In  Fig.  4,  Eq.  (16)  is  graphically  compared  with  results  using  Abramovich's 
method  (numerically  evaluated  in  Ref.  8).  The  lack  of  reliable  experimental 
data  on  this  subject  precludes  any  conclusions  as  to  the  validity  of  Eq.  (16) 
or  Abramovich's  method.  However,  in  Ref.  9,  Abramovich's  method  is 
indirectly  shown  to  be  in  general  agreement  with  experiment. 
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The  width,  b,  of  a  turbulent  mixing  zone  can  be  derived  in  the  same 
manner  as  for  the  laminar  mixing  case.  The  basic  equation  is 


b  =  (£) 

Substituting  Eq.  {15)  into  Eq.  (17)  yields 

b  _  (4f)D  * 

(I)>7d  (!  “  c.^  ^d) 


(17) 


(18) 


The  mixing  zone  width  based  on  the  usual  definition  of  the  mixing 
region,  =  1.503  ,  is 


(  0.1273 )  (4f)p  X 
(1*^D  0  "  C*o0  ) 


(19) 


Thus,  for  given  free-stream  conditions,  the  width  of  a  turbulent  mix¬ 
ing  zone  varies  linearly  with  X  .  The  q  ordinate  for  turbulent  mixing  is 


C1  “  Cg-  )  /y\ 

(0.04235)  (4|-)d  W 


SECTION  III 
CONCLUSIONS 


The  velocity  profile  for  laminar  or  turbulent  two-dimensional  mixing 
is  approximated  very  well  by  the  equation 


$  =  -1—  [1  +  erf  q  ] 


(21) 


where 


a  Y 
X 


Within  the  framework  of  the  assumptions  of  this  analysis,  the  similarity 
parameter,  <7,,  for  laminar  mixing  is 


(22) 
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The  location  of  the  dividing  streamline  for  laminar  mixing  has  been 
theoretically  determined  by  Chapman  and  is  independent  of  free-stream 
velocity.  Thus,  the  laminar  mixing  problem  is  theoretically  defined,  in 
closed  form,  by  Eqs.  (21)  and  (22)  and  Chapman's  results  without  the 
need  of  a  reference  perturbation  velocity  factor,  a  . 


The  similarity  parameter  for  turbulent  mixing  based  on  Prandtl's 
mixing  length  theory  is 


(0.04235)  (■*£)' 


(23) 


The  numerical  constant  was  determined  from  the  experimental  results 
obtained  by  Tollmien.  However,  this  constant  is  independent  of  the  defini¬ 
tion  of  the  mixing  region  and,  based  on  the  agreement  with  Abramovich's 
method,  also  independent  of  the  free-stream  velocity.  These  results  con¬ 
firm  the  applicability  of  Prandtl’s  mixing  theory  to  this  type  of  compres¬ 
sible  mixing, 


The  similarity  parameters  for  both  laminar  and  turbulent  mixing 
were  derived  from  a  single  differential  equation,  which  relates  u  to  the 
mixing  length,  X,  and  the  viscous  shear  stress  on  the  dividing  streamline, 
vD  .  This  equation  is 


O 

p» 

h— i 

■w--1 

"  8 

9 

QL 

-(*) 

tD  - 

1 

CC 

_ 1 

d  X 

(24) 


Although  this  is  an  analysis  of  a  relatively  simple  mixing  zone,  the 
approach  can  be  applied  to  more  complex  mixing  processes  such  as  two- 
dimensional  or  axisymmetric,  non-isoenergetic  mixing  with  or  without  a 
secondary  stream. 


REFERENCES 


1.  Chapman,  D.  R.  "Laminar  Mixing  of  a  Compressible  Fluid.  " 

NACA  Report  958,  1950. 

2.  Nash,  J.  F.  "Laminar  Mixing  of  a  Non-Uniform  Stream  with  a 

Fluid  at  Rest.  "  ARC-TR-22-245  or  CP-613,  September  1960. 

3.  Pai,  S.  I.  Fluid  Dynamics  of  Jets.  D.  Van  Nostrand  Co. ,  Inc., 

New  York,  1954. 


10 


AEDC-TR-65-84 


4.  Korst,  H.  H. ,  Chow,  W.  L. ,  and  Zumwalt,  G.  W.  "Research  on 

Transonic  and  Supersonic  Flow  of  a  Real  Fluid  at  Abrupt  In¬ 
creases  in  Cross  Section.  "  ME-TR-392-5,  University  of 
Illinois,  December  1959. 

5.  Page,  R.  H.  and  Dixon,  R.  J.  "A  Transformation  for  Wake  Analyses.  " 

A1AA  Journal,  Vol.  2,  p.  1464,  August  1964. 

6.  Korst,  H.  H. ,  Page,  R.  H. ,  and  Childs,  M.  E.  "Compressible  Two- 

Dimensional  Jet  Mixing  at  Constant  Pressure,  Tables  of  Auxiliary 
Functions  for  Fully  Developed  Mixing  Profiles.  "  ME-TN-392-3, 
University  of  Illinois,  April  1955. 

7.  Tollmien,  W.  "Calculation  of  Turbulent  Expansion  Processes.  " 

NACA-TM-1085. 

8.  Bauer,  R.  C.  "Characteristics  of  Axisymmetric  and  Two-Dimensional 

Isoenergetic  Jet  Mixing  Zones. 11  AEDC-TDR-63-253  (AD426116), 
December  1963. 

9.  Bauer,  R.  C.  "Theoretical  Base  Pressure  Analysis  of  Axisymmetric 

Ejectors  without  Induced  Flow."  AEDC-TDR-64-3  (AD428533), 
January  1964. 


11 


UNCLASSIFIED 
Security  Classification 


DOCUMENT  CONTROL  DATA  •  RAD 

_ (Security  clarification  of  title,  body  ot  ebmtract  and  indexing  anno  fatten  mu&t  be  entered  when  die  overall  report  i§  claesilied) 


1.  ORIGINATIN  G  ACTIVITY  (Corporate  author) 

Arnold  Engineering  Development  Center 
ARO,  Inc.  Operating  Contractor 
Arnold  AF  Station.  Tennessee 


3.  REPORT  TITLE 

AN  ANALYSIS  OF  TWO-DIMENSIONAL  LAMINAR  AND  TURBULENT  COMPRESSIBLE 
MIXING 


2a-  REPORT  SECURITY  C  CLASSIFICATION 

UNCLASSIFIED 


23  GROUP 

N/A 


4  DESCRIPTIVE  NOTES  (Type  of  raport  a nd  Inclusive  d a  lea} 

N/A 


5-  AUTHORfSJ  (Laet  name,  II ret  name,  initial) 

Bauer,  R.  C.,  ARO,  Inc. 


6-  REPORT  DATE 

April  1965 


8<l  CONTRACT  OR  GRANT  NO. 

AF  40 (600) -1000 

b.  PROJECT  NO. 

6950 

'■Program  Element  62405184 
4. Task  695002 


7a-  TOTAL  NO.  OF  PAGE*  73.  NO.  OF  REPS 

17  9 


•  a.  ORIGINATOR’*  REPORT  NUMBERfSj 

AEDC-TR-6  5-84 


S3.  OTHER  REPORT  NOfS>  (Any  other  number,  (hat  may  3a  aaaldnad 
ft  (a  report) 

N/A 


10-  A  VA  IL  ABILITY/LIMITATION  NOTICES 

Qualified  requesters  may  obtain  < 

copies  of  this  report  from  DDC . 

11.  SUPPLEMENTARY  NOTES 

N/A 

12.  SPONSORING  military  activity 

Arnold  Engineering  Development  Center 
Mr  Force  Systems  Command 
krnold  AF  Station,  Tennessee 

is.  abstract  An  analysis  is  presented  of  two-dimensional,  isoenergetic , 
compressible  mixing  of  a  jet  with  a  fluid  at  rest  for  both  laminar 
and  turbulent  mixing.  The  analysis  is  primarily  concerned  with  the 
development  of  theoretical  expressions  for  the  mixing  similarity 
parameters.  A  general  momentum  equation  is  derived,  which  relates 
the  similarity  parameter  to  the  mixing  length  and  viscous  shear 
stress.  By  using  this  general  equation  and  Newton’s  viscous  shear 
relation,  a  complete  theoretical  solution  for  the  laminar  mixing 
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regions  are  derived.  These  results  and  the  approximate  theoretical 
velocity  profile  equation  developed  by  Pai,  Nash,  and  Korst  represent 
a  closed-form  theoretical  approximation  of  the  type  of  mixing  con¬ 
sidered. 
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